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Abstract. - We investigate the trapping problem in Erdos-Renyi (ER) and Scale- Free (SF) net- 
works. We calculate the evolution of the particle density p(t) of random walkers in the presence 
of one or multiple traps with concentration c. We show using theory and simulations that in ER 
networks, while for short times pit) oc exp(— Act), for longer times p(t) exhibits a more complex 
behavior, with explicit dependence on both the number of traps and the size of the network. In SF 
networks we reveal the significant impact of the trap's location: p(t) is drastically different when 
a trap is placed on a random node compared to the case of the trap being on the node with the 
maximum connectivity. For the latter case we find p(t) oc exp for all 7 > 2, where 

7 is the exponent of the degree distribution P(k) oc k 



Introduction. — The properties of random walk greatly vary depending on the di- 
mension and the structure of the medium in which it is confined [1-4] , where a particularly 
interesting medium for the study of the random walk is complex networks [5-9] . Networks 
describe systems from various fields, such as communication (e.g. the Internet), the so- 
cial sciences, transportation, biology, and others. Many of these networks are scale-free 
(SF) [10-13]. This class of networks is defined by a broad degree distribution, such as a 
power law P{k) oc fc -7 (k > m), where 7 is a parameter which controls the broadness of the 
distribution. 

Trapping is a random walk problem in which traps are placed in random locations, 
absorbing all walkers that visit them. This problem was shown to yield different results 
over different geometries, dimensions and time regimes [2,3,14-17]. The main property of 
interest during such a process is the survival probability p(t), which denotes the probability 
that a particle survives after t steps. The problem was studied in regular lattices and in 
fractal spaces [2,14-19] and recently, in small- world networks [6]. 

In this Letter we study the problem of trapping in networks. This is a model for the 
propagation of information in certain communication networks. This follows since in some 
cases data packets traverse the network in a random fashion (for example, in wireless sensor 
networks [20], ad-hoc networks [21] and peer-to-peer networks [22]). A malfunctioning node 
in which information is lost (e.g., a router which cannot transmit data due to some failure) 
acts just like a trap in the model. This model can also be applied to loss of information 
in messages over communication systems, e.g. in the case of e-mail messages, where a 
malfunctioning e-mail server acts as a node absorbing, but not transmitting, all e-mail 
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messages it receives. Furthermore, our model may be relevant in social systems, where some 
information may initially spread randomly, but in later stages it might be held by certain 
individuals. 

We study the survival probability p(t) of random walkers on random regular networks 
(networks in which all nodes have equal degree), Erdos-Rcnyi (ER) networks (a simple 
model for random networks in which all links exist with the same probability [23-25]), and 
SF networks. We derive analytical expressions for p(t) for a wide range of trap degrees 
and concentrations and highlight the role of the network structure, obtaining new scaling 
relations for the survival probability and average trapping time which are absent in lattices. 
Our analytical predictions are confirmed with Monte-Carlo simulations. 

Methods. — To perform Monte-Carlo simulations, we generate ER networks by con- 
sidering all pairs of nodes and linking a pair with probability p. The construction of an SF 
network follows the Molloy - Reed scheme [26] . Each node i is assigned a number of links 
taken from the distribution P(k) cx fc~ 7 and then open links are connected randomly. The 
value of k is taken to be between m (typically 1-3) to k max = N — 1 (no upper cutoff value 
is imposed). We find the largest cluster by using depth- first search [27] and then discard all 
nodes that are not in the largest cluster. Starting from a fixed density of particles initially 
placed in random nodes, particles hop with equal probability to one of their nearest neigh- 
bors. Certain nodes are randomly chosen to serve as traps. These are perfect traps; if a 
particle falls on it then it is trapped and removed from the network. In the case of multiple 
traps, n = cN traps are placed in the network, where c is the trap concentration. 

Results. — Assume the network has N nodes, average degree (k) and n traps. How 
does p change as t increases to t + 1 (i.e., after each particle has moved once)? Denote 
the traps by (ii,i2, ■■■in) an d define k n = k ix + k i2 + ... + k in as the total number of links 
emanating from all traps. If at time t a given particle is not on a trap, but will hop on 
any of these k n links on its next step, it will be trapped at time t + 1. We approximate 
the probability for the particle to hop on any of these k n links to be proportional to their 
relative number in the network, that is, Ajffhr, where N (k) is the total number of links 
in the network, and A = 0(1) is the proportion constant which we will study later. In 
continuous time, this results in the equation: 

dp . k n . . 

dt = - pA Njk) (1) 

whose solution is: 

p(t)=p exp[-Atk n /(N(k))] (2) 

Surprisingly, although Eq. (2) is based on a rather simple approximation, we show below 
that it predicts very accurately the survival probability for various network models, time 
scales, and trap concentrations. In fact, Eq. (2) can be seen as a special case of the theory 
developed earlier in [15,17,28,29], where it was shown, that for a ci-dimensional lattice, 
the survival probability decays as a stretched exponential p ~ e~ at with [3 — -A^- Since 
networks have infinite dimension, d/(d+2) — > 1 to recover the exponential decay we predict. 
Note that the average time before trapping is O(N), as expected from first passage time 
considerations [30] (see also below). 

A necessary condition for the above approximation to hold is that the number of links 
between the traps is negligible. For ER networks where links exist independently of one 
another the probability that all links emanating from the traps connect to non-traps nodes 
is [(N - n)/(N - l)f". Since (k n ) = n ( k), and [1 - n/(N - l)]" <fe> ss 1 - n 2 (k) /N (for 
n <C N), we expect that as long as n -C y/N/ (k), this condition is satisfied. 

In the following, we will apply Eq. (2) to specific network topologies. In random regular 
networks, where each node has exactly k neighbors, we use Eq. (2) by substituting (k) = k 
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and k n — nk: 

p = exp [— Ant/N] = cxp [—Act] (Regular, n traps), (3) 

(without loss of generality, we set p(0) = 1). For other networks one has to take into account 
the distribution of degrees. Thus, in order to average p{t) over all networks in the ensemble, 
we need to condition p on k n : 



P(*)net = ]T p { fc « = k}exp[-Atk/(N(k))]. 



(4) 



Consider ER networks with one trap (n = 1): P{k n = k} = e W (k) /fc! is the degree 
distribution (a Poisson) [23-25]. Thus: 

=£>-<*> 
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(5) 



(ER, one trap, random k) 



where we start the summation from k = 1, since we do not place a trap on an isolated 
node (k = 0). However, when evaluating the sum, we assume the probability for k = 
is negligible, which is justified whenever (k) is large enough, which we assume henceforth. 
Also, in the simulations, we consider only the largest connected cluster, which by definition 
contains no isolated nodes. 

The same approach can be applied to the case of multiple traps, by realizing that (ne- 
glecting links between the traps) the sum of links emanating from the traps is a sum of 
Poisson variables with mean (fe), which is itself a Poisson with mean n (k): 



P(t) = jr e 



ra(fc) 



k=n 



(n(k))> 
k\ 



exp 



Atk 



■ exp 



— n (k) ( 1 — cxp 



N(k) 
At 



N(k) 

(ER, n traps, random fe) 



(6) 



The agreement of Eq. (6) with simulation results is evident from Figure 1. Note that the 
survival probability in Eq. (6) does not solely depend on the trap concentration c = n/N, but 



At 



N(k) 



At 



N{k) 



on both n and N, except for the short time limit (t <C N (fc)), when 1 — exp 

and p rts cxp [— Ant/N] = cxp [—Act]. For long times (t » N (fc)), due to the exponential 
dependence on t, the main contribution to the survival probability comes from configurations 
in which k n is small, the probability of which depends on n alone. On the other hand, the 
probability that the particle falls into the trap still depends on the total number of links 
N(k). Thus, the survival probability depends on both n and N independently. It can 
also be seen that particles survive longer as the network becomes smaller (Figure 1(b)) and 
sparser (Figure 1(c)). 

Even though scale-free networks are highly heterogeneous and thus the approximate 
approach is expected to yield less accurate results, nevertheless it is still quite useful. The 
degree distribution is P{k) = Ck~ 7 , k > to, where C is a normalization factor. Thus, for a 
single trap: 

i o = ^Cfc- 7 exp[-Atfc/(7V(fc))]. (7) 

k 
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Fig. 1: Trapping in ER networks, (a) Particle density p(t) vs. t (measured in Monte-Carlo steps). 
The network parameters are: N = 10000 and (k) = 10. Traps are placed with a concentration c on 
random nodes of the network. All results are averaged on at least 5000 runs, each with a different 
configuration of the network. Solid lines represent fitting with Eq. (6) (with the number of traps 
n = cN). (b) p(t) for fixed trap concentration c = 0.001, average degree (k) = 10, and different 
system sizes, (c) p(t) for fixed trap concentration c = 0.001, system size N = 10000, and different 
average degrees. 



Since this docs not lead to a closed form formula, we focus on the case where the degree of 
the trap k is fixed. We expect: 

p = exp [-Akt/(N (k))\ (SF, one trap, fixed k). (8) 

Interestingly, simulations show a distinct behavior for m < 3, and m > 3 (Figure 2). While 
in the case of m > 3 the simulations agree with the theory (Eq. (8)), as is evident by 
the collapse of all curves with the same kt; for m < 3 the decay of p(t) is slower than 
exponential. Note that in contrast to ER networks, p(t) is larger for the denser networks 
(smaller 7). Thus, whereas ER networks become less robust as links are added, SF networks 
gain robustness. This is a fundamental difference between ER and SF networks revealed by 
our results. 

When the degree of the trap is allowed to vary, we consider the long time regime. As in 
ER networks, the main contribution comes from configurations in which the degree of trap 
is minimal i.e., k = m. Thus we expect: 

p(t) w cxp [-Amt/(N (k))} 

(SF, one trap, random i,m>3,t>iV (k)) (9) 

which agrees with simulations (see Figure 3(a)). For SF networks with many traps, a 
simple generalization of Eq. (9) (replacing m by nm) is not applicable, and we report only 
the numerical results (Figure 3(b)). Here, similarly to ER networks, the smaller networks 
are more robust. 

SF networks exhibit nodes of particular importance which have many connections and 
play special role in transport [31]. Thus, it is interesting to study a failure in the node of 
highest degree (the hub) [32], which results in trapping of incoming particles. The maximum 
degree K in SF networks scales like K w mA^~ (for 7 > 2) [33]. Substituting k = K in 
Eq. (8), we find: 

p = cxp[- AtK/N (k)} ss exp -Amt/N^ (k) 

(SF, trap on the hub, m > 3, t > N (k)) (10) 
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Fig. 2: Trapping in scale-free networks with a single trap on a node with fixed degree k. (a) Particle 
density p(t) vs. kt, for SF networks with N = 10000, 7 = 2.5, m = 1, 2, and different trap degrees, 
(b) Same as (a), but for m = 3 (7 = 2.5, 5). In this case, all curves collapse, in agreement with Eq. 
(8). 



and the average time before trapping is thus t tI ~ N 13 where (3 = ^—j < 1 (see Figure 4). 
Realistic SF networks have 2 < 7 < 3 [10,34,35], so that < (3 < 1/2. This implies that real- 
world networks are ultra-prone to failure in their highly connected nodes. This is an even 
stronger effect compared to the targeted removal of high degree nodes [32] , whereas a failure 
of only one hub induces a significant decrease in the trapping time, a finite concentration of 
hubs has to be removed to fragment the network. 

The dramatic decrease in the time before failure is not limited to placing the trap pre- 
cisely on the node with the maximal degree. It can be proven that whenever we either (i) 
choose the node with maximal degree out of n random nodes, when n/N is finite, or, (ii) 
choose one of the nth nodes of highest degree when n = 0(1), the probability of the trap 
degree to exceed K — miV~ is finite. Thus, in these cases, the trap will be attached to a 
sufficient number of links for the scaling t tr <~ to appear. 

The value of 7 for which SF networks are equivalent to ER networks is a topic of recent 
interest [36]. Our results suggest that SF networks are equivalent to ER only when 7 
is infinite, since only when 7 — > 00 does /3 — > 1, as for homogenous ER networks. For 
ER networks, the degree distribution is a Poisson with variance equals to the mean (k). 
Consequently, the typical maximal degree is roughly K w (k) + (k). This yields p w 



exp 



At 
N 



1 + 



such that the typical time is t tI ~ as before. 
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Fig. 3: Trapping in scale-free networks with m = 3. (a) Particle density p(t) vs. t, for SF networks 
with 7 = 2.5, a single trap on a random node, and different system sizes. Solid lines represent 
fitting to exponential decay in the long time regime, Eq. (9). (b) Particle density p(t) vs. t, for SF 
networks with 7 = 2.5, traps with concentration c = 0.01 placed on random nodes, and different 
system sizes. 



In the following, we study the behavior of the prefactor A. For fully connected network 
and large N, A — > 1 (sec, e.g., [37]). For sparse networks where the particle might be far 
from the trap, A is less than one, reflecting the fact that the probability to follow a link 
to the trap is somewhat less than k n /(N (k)). To find the value of A, we first point out 
that the trapping problem is a special case of a first passage time problem [3,4,8] (since 
p(t) = 1 — Y^f where F(t) is the probability to reach the trap for the first time at 

time t). To calculate the first passage time in networks, Baronchelli and Loreto [38] used 
an approximate method that exploits the small- world nature of most networks [34,39]. In 
theory, using the adjacency matrix one can calculate the transition probability matrix of 
the random walker, from which the first passage time can be easily obtained via consecutive 
powers of the matrix (see below). However, this is not feasible for large networks, and 
thus the original random walk process was reduced to a random walk between the network 
layers [38]. Given the trap, the number of nodes ri£ that are in distance I from it is calculated. 
Then, a matrix B of size ^ max x ^max is constructed, in which Bg^i is the probability of a 
random walker in layer i to jump into a node in layer £'. For most real and model networks, 
f max ~ log N such that the size of the problem is reduced exponentially. Define the number 
of links that connect layers I and I + 1 by si, the number of links within layer t by 01 , and 
the sum of degrees of nodes in layer £ by ran = {se + se-i + 2o^). Since the random walker 
jumps into each link with equal probability, the only non-zero elements are: -B^+i = Si/me, 
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Fig. 4: Trapping in SF networks after failure of the most connected node, (a) Particle density p(t) 
vs. t, in SF networks with 7 = 2.5 and m = 3, for different system sizes. One trap is placed on 
the node with the maximum degree. Solid lines represent fitting to an exponential decay p ~ e~ Bt 
in the long-time regime, (b) The exponent B vs. N, for 7 = 2.5 plotted in (a), as well as for 
7 = 3,4, 5. It can be seen that B ~ N^ 13 with (3 « (inset), in agreement with Eq. (10). 



Bi^t-i = Si-i/me,, and Bi t i = 2o£/m,£. To represent the trap, Bo,e = for all I (since the 
trap forms layer i = 0). Since the probability for a random particle to start in a node of 
layer i is ne/(N — 1), the first passage time probability F(t) is given by [38]: 

^max 

In [38], the matrix B was constructed for ER networks, and it was found that F(t) cx 
exp [— Akt/N (fc)], with the value of A determined numerically. In the following, we extend 
this approach to study random regular and SF networks which were not previously studied, 
and the behavior of A for a concentration of traps and different average degrees in ER 
networks. 

To construct the transition probability matrix for n traps in ER networks, one has 
to calculate the number of nodes within distance I from each of the n nodes. This is 
easily accomplished by setting no = n and using the same formulae as in [38] n^+i = 

(n — J2k=o Uk ) [1 ~ (1 ~ P) ne ]i w ith P = (k) /(N — 1) being the independent probability of 
a link to exist. The number of links connecting the different layers is the same as in [38]: 
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Fig. 5: An analysis of the prefactor A. (a) Plot of 1 — A vs. (k) in ER networks, where A is the 
fitting parameter in Eq. (5). Here N = 10000 and one trap was placed randomly (black squares). 
The solid line corresponds to the theoretical result A = 1 — 1/ (k) derived in the text, (b) A vs. 7 
for SF networks with N = 10 8 and m — 3, when the trap is located at the hub (theory only). 



= ni — X)fc=o Uk j P> an< ^ ° e = n e( n e — typ/^i fr° m which the matrix B is determined. 
We then calculate F(t) from Eq. (11), and A by fitting F(t) to an exponential. We find 
that the change in A for n > 1 is minor (of the order of 0(N~ 1 j) and proportional to n: 
A(ri) — A{n = 1) oc n. 

We next study the dependence of A on the average degree (for a single trap). Applying 
the above method for ER networks with different average degrees, we find that Am 1 — 1/ (k) . 
This is also confirmed by simulations (Figure 5(a)). For random regular networks, we derived 
the transition matrix B, from which we found that A w 1 — l/(fc— 1). Curiously, both results 
can be written as A w 1 — 1/(k — 1), where k — 1 = (A: 2 ) / (k) — 1 is the branching factor 
of the network (since in ER networks k = (k) + 1, and in regular networks k — k). A 
qualitative explanation of this relation (which can also be recast as A w 1 — p c , where p c is 
the percolation threshold [33]) is still lacking. 

For SF networks, we calculate A for the case when the trap is located at the node of 
maximal degree, by using [40] for the number of nodes in layer £ (ni), and the number of 
links emanating from layer I into itself (pi) and into layer £+1 (sg). As before, construction 
of the transition matrix B, application of Eq. (11), and fitting to an exponential are used to 
calculate A for different values of 7 (Figure 5(b)). Using this method we can predict A for 
very large A^s in which simulations are not possible. Here the relation j4«1-1/(k-1) [33] 
is valid up to 7 < 4. 

Conclusions. — We study the trapping problem on regular, ER, and SF networks 
using theory and simulations. We develop a simple theory to account for the behavior of 
the survival probability in a variety of conditions. In ER networks we find that the trapping 
process exhibits a non-exponential behavior which depends on both the number of traps 
and the size of the network. For SF networks we find anomalous behavior for networks with 
small minimal degree, expressed as deviations from the theory. We also find that as opposed 
to ER networks, particles survive for longer times in denser SF networks. Finally, when the 
trap is placed in one the network hubs, we find a new scaling with the system size. The 
average time before trapping decreases dramatically in comparison to random failure or to 
ER networks. This is true for all values of 7, suggesting that the equivalence of SF and ER 
networks for 7 > 4 does not exist for the trapping problem. 

Appendix 1: Finite probability for a choice of high degree traps. We prove 
the following. Given a scale-free network with N — > 00 nodes and degree distribution 
P(k) = m^- 1 ( 7 - l)fc-T , k>m: 

Theorem. 

1. The probability that the node with maximal degree out of n = cN random nodes will 



p-8 



Trapping in complex networks 



have degree that exceeds mN-i- 1 is finite, provided that c = n/N is finite. 

2. The probability that the degree of the node of nth largest degree will exceed mN^ 1 
is finite, provided that n = 0(1). 

Proof. 

1. The probability of the node with maximal degree out of n random nodes to have 

r u i cN 

degree K is V(K) « P(K) Q P(k')dk' . Substituting P(k) = cN(<y- l)™?- 1 ^ 
we have: 

V(K) = cN(~f - \)m?- x K-^ [1 - (K/m) 1 ^] cN . (12) 
The probability for K to be at least mN~ is: 



P{K > mN~} = 



^_ V(K)dK 

mN-i- 1 



roo 

/ ciV( 7 - ^trI^K-i [1 - (K/m) 1 ^} dK 

JmNi- 1 



> cN(j- l)m'>- 1 I _^K^ 



1 — (mN-i- 1 jrn) 



l- 7 



dK 



poo 

cN(j - l)™ 1 " 1 / K^e- c dK 
dV( 7 - \)rn<- x e- c [Jf 1 -T/(1 - 7)] °° 



■mNi- 



= ce 



(13) 



Thus, P-fif > mN-1- 1 } > ce c . Since c is finite, in every choice of cN nodes there is 
a finite probability that at least one of them will have degree larger than mNi^ 1 . 

2. The probability of the node with nth largest degree in the network to have degree 



K is V(K) = ( N n )nP(K) [/~ P(k')dk'] n ~ 1 [£P{k')dk> 
distribution we have: 



N-n 



. Inserting the degree 
(14) 



1-7I N ~ n 



V(K) = (^j n(7 - \)Tni- x K-~i [(K/m) 1 ^] " 1 [l - (K/r 

1 

The probability for K to be at least mN-r- 1 is: 

/>OC 

P{K >mNT^} = / t V(K)dK 

JmNi- 1 

t\ roc 

)n / _^( 1 -l)m' 1 - 1 K-"<[(K/m) 1 -" / ] n ~ 1 [l-iK/mf-^^dK 

'■/ JmNi- 1 



> 



nm "(7-l)( 7 _ 1) 



|£-n(l-7)-l 



mNi- 



= f Ar )nm"^- 1 )( 7 -l) f t A^(i-7)-i e -i 



1 - (toA^— /m) 1 " 7 



JV 



nm 



"(7-1) ( 7 _ l) e -l _ 7 )]' 



J K=mNl~ 1 



(15) 



where in the last step we used Stirling's approximation. Thus, if n = 0(1), the probability 
that the node with nth largest degree is greater than mNi^ 1 is finite. 
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Appendix 2: Distances distribution in random regular networks. We derive 
the following results for the distribution of distances from a random node in random regular 
networks. We imagine a process in which all nodes have initially k open links, and as the 
algorithm proceeds, we connect open links to form the network edges. At each step, we 
connect the open links emanating from layer £ into randomly open links from nodes which 



are outside layers 1, ...£ to form layer £+1 (see [40 
£ from the root as rig. Define rg = 1 - 



). Denote the number of nodes in distance 
12e>=o ne ' 1^ as the fraction of nodes outside layers 
1, ...,£. Define also se as the number of edges connecting nodes in layer £ to nodes in layer 
£+1, og as the number of edges connecting nodes within layer £, and eg as the number of 
edges emanating from layer £ except for the edges incoming from layer £ — 1 . The following 
recursion relations hold: 



ne+i = Nr e 



ei — 1 

III'- 

3=0 

re+i = r e - n e+1 /N 



e t - j + Nkr e 



kne + i — se 

„2 



with m = 1, r\ 



ee+i 
og+i 

1 - 1/N, ei = k, si 



-£+1 



2(kNr e+1 +e l+l ) 



(16) 

(17) 
(18) 

(19) 
(20) 



k. 
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